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Abstract 

The effect of the anti-rotating terms on the short-time evolution and the quantum 
Zeno (QZE) and anti-Zeno (AQZE) effects is studied for a two-level system coupled 
to a bosonic environment. A unitary transformation and perturbation theory are used 
to obtain the electron self-energy, energy shift and the enhanced QZE or the AQZE, 
simultaneously. The calculated Zeno time depends on the atomic transition frequency 
sensitively. When the atomic transition frequency is smaller than the central frequency 
of the spectrum of boson environment, the Zeno time is prolonged and the anti-rotating 
terms enhance the QZE; when it is larger than that the Zeno time is reduced and the 
anti-rotating terms enhance the AQZE. 
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The quantum Zeno effect (QZE) and anti-Zeno effects (AQZE) have been widely 
discussed for decades theoretically [1] and recently experimentally [2] . In an unstable 
quantum two-level (multi-level) system, frequently measurement can reduce or accel- 
erate the decay processes [1, 3, 4, 5, 6, 7]. The survival probability at an excited state 
of a quantum system interacting with an environment is a decaying function of time. 
Frequent measurements at extremely short time interval may slow down the decay 
process, because the decay of the excited state is almost zero at the beginning of the 
decay process[l, 3], which is known as the QZE. It was also found that if the mea- 
surement time interval is short, but not extremely short, the decay of the excited state 
could be accelerated [3], which is known as the AQZE. Let P(r) denote the survival 
probability (after a short time interval r) at the initial state, which can be written as 
P(t) = exp(— 7(t)t). After iV time measurements at equal r, the survival probability 
reads P N (r) = exp[— 7(r)iW] = exp[— j(r)t] with 7(7-) the effective decay rate. If 
N — 1, P(t) = exp(— j(t)t), which goes to P(t) — > exp(— j t) for large enough t, where 
7o is the decay rate under the Weisskopf-Wigner approximation. We will have the QZE 
if 7(t) < 7o, and the AQZE if 7(7-) > 70. 

It is well known that the whole spectrum of the environment at off-resonance (with 
the transition frequency) is important for the QZE and AQZE. If the atomic transition 
frequency is located not at the maximum of the spectrum, we can have the QZE and 
AQZE depending on the measurement time interval. In the previous studies on QZE 
and AQZE [1, 3, 4, 5, 6, 7], the rotating wave approximation (RWA) is used. However, 
this approximation raises the question what the influence of the anti-rotating terms 
on the QZE and AQZE is, as they might have the same order contribution as the 
spectrum components off-resonant with the atomic transition, especially for the QZE 
where measurement time interval is extremely short [1, 3, 4]. Therefore we ask ourselves, 
"what is the role of anti-rotating terms on the QZE and AQZE?" 

The model to describe an unstable quantum system is the following spin-boson 
model with the Hamiltonian[8, 9], 
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H = -u (r z + Y / UJ kb\b k + -Y.dkibl + h)a x , (1) 

k k 
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b\ (bk) is the creation (annihilation) operator of boson mode with frequency Uk, <J X 
and a z are Pauli matrices describing the two-level system, ujq is the transition fre- 
quency between the up |+) and down state |— ): <x z |±) = ±|±). gu is the coupling 
between the two-level system and the environment, which can be characterized by the 
interacting spectrum[3, 4, 8]: G(u) = \ j2k9k$( UJ ~ ^fc)- This model Hamiltonian is 
used for a large number of different physical and chemical processes, such as the atom- 
field interaction and the QZE in quantum optics[3, 4, 9], coupled quantum dots on 
a solid state substrate[8, 10], and the macroscopic quantum coherence experiment in 
SQUID's[ll, 12]. 

The Hamiltonian (1) cannot be solved exactly and usually the RWA is used [3, 4, 9] 
for which the Hamiltonian is taken to be 

1 , 1 

H rwa = -u Vz + Y, UJ kblb k + - ^g k {b\(r- + b k cr + ), (2) 

k k 

where a± = \{p x ± i> a y)- Hrwa can be solved in the so-called one-boson sector with 

the initial state | -0(0) ) = |+)|{0jt}), where |{0fc}) is the vacuum state for every k. The 

survival amplitude of finding the system still in | -0(0) ) at r > is [3, 4] 

1 r €? T dp 
xrwa(t) = — -5-, (3) 

2m jB P + ^o + \Ek^ k 

where B is the so-called Bromwich path. The survival probability in the initial state 
is Prwa(j) = \xrwa(t)\ 2 and the effective decay rate ^RWAij) for a short interval r 
can be calculated as follows [3], 

POO 

1rwa(t) = 2tt / duG(uj)F(uj - u ), (4) 
Jo 



F(uj — ujq) = 2 sin 2 



UJ — UJq 

T 



/ht(uj - u ) 2 . (5) 
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Since F(u> — uo) ~^ ^(^ — ^o) (the Dirac 5— function) for large enough t, we have the 
decay rate 70 = 2irG(u ) in the Weisskopf-Wigner approximation. Eqs.(3) and (4) are 
main results of the RWA. 

When the anti-rotating terms are included, the above method is no longer valid. 
Here we present an analytical approach, based on unitary transformation and per- 
turbation theory to calculate the survival amplitude and the effective decay rate for 
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Hamiltonian (1) in order to clarify the impact of the anti-rotating terms on the short 
time evolution and on the QZE and AQZE. In the following the interacting spectrum 
of the environment is assumed as, 

G[uj) = FifTiv' (6) 

where coupling strength, a, is a dimensionless constant and Q is the center of gravity 
of the spectrum[ll, 12], and T is the width of the spectrum. We will show that the off- 
resonance ratio ujq/VL plays an important role and the effect of the anti-rotating terms 
must be taken into account especially for the off- resonant case lu /Q <^ (1 — F/Q). We 
note that when F 2 <^ Q 2 , G(uj) is mainly a sharp Lorentzian-type peak similar to the 
case of resonant Rabi oscillation [3, 11, 12]. Throughout this paper we set h — 1. 

We treat the anti-rotating terms by a unitary transformation [13]: H' = exp(S)H exp(— S) 
with 

S = T,£r*M-b k )v x . (7) 

k Uk 

Here we introduce in S a /c-dependent function and its form will be determined later. 
The transformation can be carried out, and the result is H' = H' Q + H[ + H' 2) 

h' q = \j)^o z + 5> fc &J&* - -r-W* - 60, (8) 

2 k k AuJ k 

H i = oY,9k{l - + b k )a x - -V^yJ2—^(bl-b k ), (9) 

1 k 1 k Uk 

H' 2 = \u) (T z ( cosh{^ —ik{b\ -h)}- 77 j 

2 \ k u k ) 

-\u ia y (sinh{£ ~ h)} -vY.-Ubl ~ b k )) (10) 

z \ k u k k u k ) 

where 

k k 

H' is the unperturbed part of H' and, obviously, it can be solved exactly because in 
which the two-level system and the bosons are decoupled. The eigenstate of H' Q is a 
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direct product: |±)| {rifc}), where ({n^}) is the eigenstate of bosons with n k bosons for 
mode k. In particular, the ground state of H' is \g ) = |— )|{0 fe }). 

H[ and H' 2 depend on g k and are small, which are treated as perturbation. Because 
of the definition of 77 in Eq.(ll), H' 2 contains the terms of two-boson and multi-boson 
non-diagonal transitions and its contribution to physical results is {g k ) 2 and higher. 
So, if 2 can he omitted and we approximate H' H' q + H[. H[ contains the terms of 
single-boson transition and we chose as 



so that H[ is of the form 



& = (12) 



H[ = r]UjQ ~ik \b\a- + b k a + ] , (13) 

k U k 1 J 



It is easily to check that H[\go) = 0. We note that the transformed Hamiltonian H' is 
of a form similar to Hrwa but ojq and g k /2 in (2) are replaced by r/uo and r]Uog k ^ k /u k , 
respectively. 

As H[\g ) = 0, we have H'\g ) (H' q + if() |^ ) = E g \g ) with the ground state 

2 

energy = —\t]uoq — £^^fc(2 — which is lower than the ground state energy 
£RWA _ —I^ i n the RWA . It can be seen that the third term in H' (or the second 
term in E g ) does not depend on the atomic transition frequency; it is the self-energy 
of the free electron due to the vacuum fluctuations. In previous treatment [9] the self- 
energy usually needs to be calculated separately. 

The transition frequency (originally lvq) is modified to oo a = rjuo with the modifi- 
cation factor 7] due to the anti-rotating terms. Figure 1 shows r\ as a function of a 
for various values of uj . The horizontal dotted line is rj = 1 in the RWA. We note 
that (77 — l)cu is an energy shift of the two-level system resulted from the anti-rotating 
terms. By using this unitary transformation Eq.(7), we can obtain the self-energy and 
the energy shift simultaneously. 

Since \g ) is the ground state of H', the ground state of the original H is exp(— S) |<?o) • 
Then, the survival amplitude of finding the system in the initial state is x(t) = 
(+|({0fc}| exp(— iH'T)\+)\{0 k }). Since H' is of a form similar to H RW a, x(t) can be 
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calculated in the so-called one-boson sector in the same way as xrwa(t) was, 



where \4 = r]u gkC,k/^k- The survival probability in the initial state is P{r) = |a;(r)| 2 
and the effective decay rate 7(7-) for a short interval r is 



where G'{uS) = AG{ui){t]uiq) 2 / \uj + i]uj ) 2 = G{uj)f{uj). The spectrum is modulated by 
the factor f{uj) = {2u a ) 2 / {uj+u) a ) 2 . The physics of this factor is clear. It is proportional 
to 1/ (u+u a ) 2 because it comes from the anti-rotating terms. It is equal to 1 for uj — u a , 
because the decay rate at large enough time is proportional to G(cu a ). Please note that 
f(u) > 1 (or < 1) and G'(u) >(or <)G(u) for 00 <(or >)uj a . 

For the short-time limit the survival probability is quadratic in r: P(r) = 1 — r 2 /r| 
for r — * 0, which is explicitly different from the exponential decay law[3, 4, 6]. The 
quantity t z is referred to as "Zeno time" [4, 6] and can be calculated by using Eq.(15), 



while the Zeno time in the RWA is approximately t§ wa = (J °° cLujG(uj))~ 1 I 2 . One can 
check that r z lwA appears to be independent of c^o because the integrand G(u) does 
not depend on c^o- This should not be physically correct since the short-time evolution 
P(t) = 1 — t 2 /t| should depend on where the transition frequency ujq is located in 
the interacting spectrum. Our tz without RWA does depend on u>o, since G'(cu) is 
a function of u> . When the atomic transition frequency is smaller than the central 
frequency of the spectrum the Zeno time is prolonged; when it is larger the Zeno time 
is reduced, as shown in Fig. 2. Also the energy shift is dependent on the location of uoq 
in the spectrum. The smaller the ratio uq/VL or the larger the a (stronger interaction), 
the larger the energy shift will be (see Fig. 1). 

In Fig. 3, 7(7") is plotted for ujq located in the low frequency part of the spectrum, 
uj = 0.2Q. The dashed line is the result of RWA, and one can see that for extremely 




(14) 



7(7-) = 2tt / du)G'{u))F{u - 77^0), 



(15) 




(16) 
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short time (7 r < 0.01) RWA predicts the QZE but for a short time (70T > 0.01) there 
is a possibility of AQZE. However, by taking into account the anti-rotating terms, we 
only have the QZE and no AQZE. This is a general conclusion of our calculation for the 
off-resonant spectrum with uj /Q <C (1 — T/Q), which is different from that of Ref.[3]. 

Kofman and Kurizki[3] concluded that when uj is significantly detuned from the 
maximum of G(u) at Q (so that G(ivq) <C G(Q)), the effective decay rate 7(7") grows 
with decreasing r and leads to the AQZE of decay acceleration by frequent measure- 
ments. The reason can be understood by checking G(uj) and G'(uj) in Eqs.(4) and (15). 
The modification factor f(u>) due to the anti-rotating terms on the spectrum G(uj) is 
f{uj) < 1 for uj > uj a . When uj a is much smaller (or smaller) than the spectrum cen- 
ter frequency Q, the spectrum G'iyj) is greatly flattened (or flattened) compared with 
G(u), see the inset of Fig. 3 (or Fig. 4). The dephasing function[3] F{uj — uj ) is mainly 
a single-peak function with peak at uj and width ~ 1/r. Since the integrand in Eq.(4) 
is G(uj)F(uj — ouo), when ujq is far below the maximum of G(uj), one can check that 
1rwa{t) grows with decreasing r (AQZE) because F(uj — ujq) is then probing more of 
the rising part of G(u). Our result is different from Ref. [3] because the integrand in 
Eq.(15) is G'(uj)F(uj — uj a ) and 7(7-) decreases with decreasing r (QZE) since F(uj — uj a ) 
already covers the main part of G'{uj). 

In Fig.4 we plot 7(7-) for uj = 0.511 Here RWA predicts anti-QZE for 7 r > 0.09, 
but our result predicts very weak anti-QZE for the region 0.24 < 7 r < 0.58, because 
the factor f{uj) coming from the anti-rotating terms slightly flattens the spectrum (see 
the inset of Fig.4). 

For the resonant case, uo a ~ Q, our calculations with anti-rotating terms are nearly 
the same as those of RWA, because the factor f(ou) come from the anti-rotating terms 
is almost equal to 1 around the frequency uj = Q ~ uj a , and changes the spectrum very 
little. This can be understood easily since the RWA is a good approximation for the 
resonant case. 

Furthermore, for u larger than Q, we have f(cu) > 1 for uj < uj , that is to say, 
the main part of the spectrum enhanced by the anti-rotating terms. In Fig. 5 we plot 
7(7-) for uj a = 1.5f2, the region for AQZE is wider than predicted by RWA, because the 



7 



anti-rotating terms raise the peak of G'(u); see the inset of Fig. 5. 

In summary: The impact of the anti-rotating terms on the short-time evolution 
and the quantum Zeno and anti-Zeno effects is studied for a two-level system coupled 
to a bosonic environment. We present an analytical approach, based on a unitary 
transformation and a perturbation method. With this method, we can simultaneously 
obtain the electron self-energy, energy shift and the enhancement of the quantum Zeno 
or the anti-Zeno effects. The effective decay rate is calculated. The Zeno time depends 
on the atomic transition frequency sensitively. When the atomic transition frequency 
is smaller than the central frequency of the spectrum the Zeno time is prolonged and 
the anti-rotating terms enhance the QZE; when it is larger than that the Zeno time is 
reduced and the anti-rotating terms enhance the AQZE. 
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Figure Captions 



Fig.l The renormalization factor r] as a function of the coupling a. 

Fig. 2 The Zeno time as a function of the transition frequency c^o- 

Fig. 3 Decay rate 7(7") for a = 0.02, c^o = 0.2Q and T = OAQ. The dashed line is 
the result of RWA. Inset: G{uj) (dotted line) and G'{uj) (solid line). 

Fig. 4 Decay rate 7(7") for the below resonant spectrum with a = 0.02, ujq = 0.5Q 
and T = 0.4f2. The dashed lines is the result of RWA. Inset: G{uS) (dotted line) and 
G'{uj) (solid line) 

Fig. 5 Decay rate 7(7-) for the above resonant spectrum with a = 0.02, c^o = 15Q 
and T = 0.4f2. The dashed line is the result of RWA. Inset: G(u) (dotted line) and 
G'{uj) (solid line) 
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